Mon. Not. R. Astron. Soc. 000, ITHTT1 ( 20121 Printed 8 April 2013 (MN WF$L style file v2.2) 



X-ray diagnostics of chemical composition of the accretion 
disk and donor star in ultra-compact X-ray binaries 
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ABSTRACT 

Non-solar composition of the donor star in ultra-compact X-ray binaries may have a 
pronounced effect on the fluorescent lines appearing in their spectra due to reprocessing 
of primary radiation by the accretion disk and the white dwarf surface. We show 
that the most dramatic and easily observable consequence of the anomalous C/O 
abundance, is the significant, by more than an order of magnitude, attenuation of the 
K„ line of iron. It is caused by screening of the presence of iron by oxygen - in the C/O 
dominated material the main interaction process for a E « 7 keV photon is absorption 
by oxygen rather than by iron, contrary to the solar composition case. Ionization of 
oxygen at high mass accretion rates adds a luminosity dependence to this behavior - 
the iron line is significantly suppressed only at low luminosity, log(Lx) ^ 37 — 37.5, 
and should recover its nominal strength at higher luminosity. The increase of the EW 
of the K a lines of carbon and oxygen, on the other hand, saturates at rather moderate 
values. Screening by He is less important, due to its low ionization threshold and 
because in the accretion disk it is mostly ionized. Consequently, in the case of the 
He-rich donor, the iron line strength remains close to its nominal value, determined 
by the iron abundance in the accretion disk. This opens the possibility of constraining 
the nature of donor stars in UCXBs by means of X-ray spectroscopy with moderate 
energy resolution. 
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1 INTRODUCTION 

Ultra compact X-ray binaries (UCXBs) are a sub-group of 
X-ray binaries with orbital periods of less than 1 hour. 
Their small orbital periods do n ot allow for a hydrogen 
rich, main sequence donor (e.g . iRappaport fc Josa 1 19841 : 



iNelson. Rappaport. fe JossHl986T ). T he most likely scenarios 



for their formation (for details see e. g. Tutukov fc Yungelsonl 



1 19931 : llben. Tutukov. fc Yungelsonl 1 1 19951 : INelemans! I2005T ) 



predict that the donor star in such systems is a white 
dwarf (WD) or a helium star. Driven by the loss of the 
orbital angular momentum due to gravitational wave ra- 
diation, UCXBs are often observed as persistent and rel- 
atively luminous X-ray s ources with luminosities in the 
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iBildsten fc Delove!l2004h 
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Given the nature of the donor star, the accreting mate- 
rial in UCXBs should have a chemical composition consis- 
tent with the ashes of H burning (mostly He and 14 N), He 
burning (mostly C/O) or carbon burning (mostly O/Ne). 
Depending on the binary's formation channel, it can vary 
from C/O-rich to He-rich. Indeed, optical observations of 
several UCXBs, for example, 4U 0614+091, 4U 1543-624 
and 2S 0918-549 suggest accretion of C/O-rich materia l 
(INelemans et alj |2004 iNelemans. Jonker, fc Steeghj |2003 : 
I Werner et alj 120061 b On the other hand, in the case of 4U 
1916-05 they reveal evidence poin ting to a He-rich donor 
|Nelemans. Jonker fc Steeghj|20o9 ). Modeling of type I X- 
ray bursts from 4U 1820-30 suggest that the accreting mate- 
rial in this syste m is also helium dominated (|Bildstenlll995l : 
ICumming|[2003h . 

X-ray spectra of X-ray bina ries usually con tain the so 
called reflected component (e.g. iGilfanovl |2010| . and refer- 
ences therein). This component is produced due to repro- 
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cessing of primary emission by the optically thick Shakura- 
Sunyaev accretion disk and by the surface of the donor star 
facing the compact object. The primary emission may origi- 
nate in a hot optically thin corona, in the accretion disk itself 
or, in the case of a NS accretor, in the boundary layer on the 
surface of the star, and carries most of the energy. Depend- 
ing on its origin, the spectrum of the primary emission may 
vary from soft thermal to hard power law-like spectrum. Al- 
though the reflected component is energetically insignificant, 
it carries information about the geometry of the accretion 
flow (e.g. iGilfanov. Churazov. fe Revnivtsevlll999h and, via 
fluorescent lines and absorption edges of metals, about the 
chemical composition, ionization state and kinematics of the 
accretion disk material. 

Reprocessing of X-ray radiation by the accretion disk 
and by the surface of the donor star has been stud- 
ied extensively by many authors, starting fr o m th e sem- 
inal paper by iBasko. Sunvaev" fc Titarchukl || 19741 ). The 
shape and strength of the iron K a fl uore s cent l ine has 
been investigated bv iBaskd (| 19781 ) and Ha] (jl979l ). Semi- 
analytical expressions for reflection spectra including the 
effects of both photoioni zation and Compt o n sc atter- 
ing, have been derived bv iLightman fc White! (| 19881 ) and 
IWhite. Lightman. fc Zdziarskil (|l988f ). and formulated in 
terms of K-shell fluorescence and the characteristic Compton 
hump between 10 keV and 300 keV. Later on, the effect of 
ionization of the ac creting material on reflected spectra has 
been included (e.g. IRoss fc Fabianlll993l ; IZvcki et al.lfl99l ; 
iNavakshin. Kazanas. fc K allman 2000). The need for better 
accuracy and more realistic and complex geometries has led 
to application of Monte-Carlo (MC) methods that compli- 
mented and enhanced analytical calculations. Many authors 
have comput ed detailed models based on Monte-Carlo tech- 
niques (e.g. Ge orge & Fabian 199 l]; Matt. Fabian, fc Ross! 



Table 1. Abundances of elements for different types of white 
dwarfs used throughout the paper. 



ll993l ; lBallantvne. Ross, fc Fabianll200ll ) 

Despite of the amount of effort invested in studying re- 
flection of X-ray emission from optically thick media, all 
prior work concentrated on the ~ solar abundance case, 
with only moderate variations of the element abundances 
considered in some of the papers. On the other hand, in the 
case of UCXBs, we expect that accreting material may have 
significantly non-solar abundances, for example with all hy- 
drogen and helium being converted to carbon and oxygen. 
Such drastic abundance modifications should lead to strong 
changes in the properties of the reflected spectrum, espe- 
cially in its fluorescent line content. This problem is investi- 
gated in the present paper. The composition of the accreting 
material is discussed in Section 2. In Sections 3-4, we con- 
sider an idealized case of an optically thick slab of neutral 
material in order to identify main trends and then (Section 
5) discuss modifications to this picture which may be in- 
troduced by gravitational settling of heavy elements in the 
white dwarf envelope and ionization of the accretion disk 
material by viscous heating and irradiation. We use simple 
analytical calculations (Section 3 and Appendix) to illus- 
trate the physical origin of the main dependencies and then 
utilize the Monte-Carlo technique to compute reflected spec- 
tra for strongly non-solar abundances of the type expected 
in the donor stars in UCXBs (Section 4). We mainly focus 



Element 


He 


c/o 


O/Ne 




Pandey 


Garcia Berro 


Gil-Pons &, 




et.al. 


et.al. 


Garcia Berro 


H 


1.99 ■ 10~ 6 


- 


- 


He 


0.997 


- 


- 


C 


1.58-10~ 3 


0.563 


- 





3.9710~ y < 


0.422 


0.649 


Ne 


3.66-1CT 4 


1.37-10- 2 


0.262 


Na 


6.06-1CT 6 


2.08 


4.93-10~ 2 


Mg 


1.08-HT 4 


3.69 1(T 4 


3.88-10~ 2 


Si 


1.0MCT 4 


3.44-10~ 4 


4.43-10" 4 


s 


4.59-1CT 5 


1.57-10- 4 


2.02 10" 4 


Ar 


1.27-1CT 5 


4.34-10~ 5 


5.5840~ 5 


Cr 


1.33-1CT 6 


4.54-10~ 6 


5.8440~ 6 


Mn 


6.94-1CT 7 


2.38-10" 6 


3.06 lO" 6 


Fe 


9.18-lCr 5 


3.14-10- 4 


4.0440" 4 


Co 


2.36-1CT 7 


8.07-10- 7 


1.04 


Ni 


5.04-1CT 6 


1.73-10- 5 


2.22 10~ 5 


Cu 


4.59-lCT 8 


1.57-10" 7 


2.02 10" 7 


Zn 


1.1340" 7 


3.86-10" 7 


4.96 10" 7 



Abundances are by number of particles. The references to the 
original abundance calculations are given in the column titles. 
The numbers written in italics are from these calculations. Abun- 
dances of other elements were fixed at the solar values in mass 
units and then converted to concentration abundances. 

on the strengths of fluorescent lines of the elements expected 
to be abundant in different formation scenarios of UCXBs. 



2 COMPOSITION OF THE ACCRETING 
MATERIAL 

Different initial parameters and the environment 
of UCXB progenit ors may lead to a variety of 
donors (e.j 



1986: 



Savoniie. de Kool. fc van den Heuvell 

iPodsiadlowski. Rappaport. fc Pfah] |2002| ; 



Yungelson. Nelemans. fc van den Heuvell 200^ ) 



non- 
degenerate He star, He white dwarf, C-O or O-Ne-Mg white 
dwarf. As it will become clear later, from the point of 
view of classification of the reflected spectra, the variety of 
abundance patterns can be broadly divided into two types 
- (i) He-rich and (ii) C/O/Ne/Mg-rich. 

A typical isolated C/O white dwarf is expected to con- 
sist of a core mostly made of a mixture of carbon and oxygen, 
surrounded by a He-rich layer of up to 10~ 2 Mq. On top of 
the helium layer there may be a thin layer of hydrog en of up 

A late 



tne helium layer there ma y be a tnm layer 01 Hydros 
to 10~ 4 M (e.g. lKawalerlll995l ; lAlthaus et alJhoiof ) 



shell flash could produce a C/O-rich envel ope compri sed of 
~ 30% of He from convective shell burning (Ibcn 1983). The 
mass of the He-C-O mantle surrounding the C/O core can 
grow up to 10 _1 Mq, itself envelope d by a thin layer of H 
in th e case of a hybrid white dwarf |lben fc Tutukovlll985l . 
Il987t ). On the other hand, if the initial mass of the compan- 
ion star was in the ~ 8 — IIMq range, a UCXB with an O-Ne 



white dwarf donor may be formed (|Gil-Pons fc Garcia-Berrol 
l200ll) . 
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Figure 1. The probability Pi(E) (cq. JTJ) for an incident photon to be absorbed due to K-shell ionization of carbon (blue dashed line), 
oxygen (black solid line) and iron (red dash-dotted line), in the solar abundance case (left panel) and for the composition of a C/O white 
dwarf as tabulated in Table [l] (right panel) . 



In the case of a white dwarf in a binary system, 
this basic structure will be modified by a co-evolution 
with the companion star. H and He layers can either be 
strip ped away during the initial stages o f binary interac- 
tion ([Kaplan. Bildsten. fc Steinfadtl 120121 ) or be gradually 
depleted due to accretion. Indeed, typical luminosities of 



10 erg/s range, implying the 
- TO" " 



UCXBs are in the ~ 10 d 

mass accretion rate in the ~ 10 -10 — 10 -8 Mo/yr range. 
At this rate, a surface layer of < 10 -2 Mq will be depleted 
within < 1 — 100 Myrs, which is (much) shorter than the 
expected life times of such systems. 

If the white dwarf donor has been completely stripped 
of its H and He layers, the chemical composition of the 
accreting material will be determined by its core. For the 
purpose of this calculation we will ignore the complexity 
of the possible abundance patterns and assume the follow- 
ing mass fractions : C — O ~ 0.49 and Ne ~ 0.02 (e.g. 
iGarcfa-Berro et al.l 120081 ) . More massive white dwarfs are 
expe cted to have cores composed mainly o f oxygen and neon 
(e.g. iRitossa. Garcia-Berro. fc Ibenlll996l ). in which case we 
assume the following composition: Ne ~ 0.28, O ~ 0.55, 
Mg ~ 0.05, Na ~ 0.06 jGil-Pons fc Garcia-Berroll200ll) . 

The He-star and evolved secondary donor scenarios 
yield an accreting material that is He-dominated, perhaps 
with some He burning products and traces of H depending 
on the phase at which th e evolved seco n dary star started 
its Roche lobe overflow. IPandev et al.l (|200lh have per- 
formed detailed calculations for the abundances of cool He- 
stars and obtained the following values: He=0.99, C=0.0052, 
N=0.0016, 0=0.0018. For an ev olved main sequence star 
companion iNelemans et al.l (|20ld ) used the Eggleton stellar 
evolution code TWIN to model the abundance pattern at the 
donor's surface. They predict material that is He-rich with 
less than 0.01 O and N and traces of C. 

For the purpose of this paper we will assume that the 
mass fractions of all other elements, not included in the 
above calculations, are equal to their solar values. This is 
equivalent to the assumption that the total masses of these 



elements in the star did not change in the course of its evo- 
lution towards a white dwarf. The corresponding mass frac- 
tions were then converted to the concentration abundances 
taking into account changed particle concentrations of H, 
He, C, O etc. Solar abundances fo r elements with Z=l-30 
were adopted from iFeldmanl dl992l). elements not l i sted i n 
this tabulation were taken from lGrevesse fc Sauvall (|l99ct ). 
The abundance patterns for different types of white dwarfs 
used throughout the paper are summarized in Table [1] In 
the Table, the values obtained in the white dwarf abundance 
calculations described above are given in italic, abundances 
of elements which mass fractions were fixed at their solar 
value are shown in roman font. 



3 QUALITATIVE PICTURE 

In this section and in section 4, we consider an idealized 
case of reflection from cold and neutral material. Although, 
at a sufficiently low mass accretion rate, this is a reasonable 
approximation for carbon and oxygen, it breaks down for 
helium in the accretion disk at any luminosity relevant to 
UCXBs. In the case of a cold white dwarf, it may hold on 
its surface. The effect of ionization is discussed in detail in 
Section 5. 

Although the strengths of fluorescent lines depend on 
the number of parameters, such as spectral index of the pri- 
mary radiation, inciden c e angle, ionization state etc . . (e.g . 
iGeorge fc Fabianl Il99ll : iBallantvne, Fabian, fc Ross! 120021 ) 
variation of these parameters within their plausible ranges 
leads to rather moderate changes in the equivalent widths. 
Much more dramatically the EWs of fluorescent lines are 
affected by changes in the relative abundances of different 
elements. As it turns out, an increase in the abundances 
of carbon and oxygen, not only changes the strengths of 
fluorescent lines of these elements, but has a much easier 
detectable effect on the iron line. Indeed, as it will be shown 
below, for the typical X-ray spectra of UCXBs, the equiva- 
lent widths of carbon and oxygen lines increase by ~ 2 — 10 
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Figure 2. The probability of absorption by oxygen and iron 
calculated at the respective K-edges, versus oxygen abundance 
Aq. The latter is expressed in solar units. 



times with respect to their solar abundance values, but still 
remain in the hardly detectable sub-eV range. On the other 
hand, the iron line drops more than 10-fold, from easily ob- 
servable ~ 100 eV to the ~ few eV range. 

The effect of elemental abundances on the strength of 
emission lines, can be illustrated by the following simple 
calculation. Let's consider an optically thin layer of mate- 
rial illuminated by a photon beam traveling along the axis 
normal to its surface. The probability that an incident pho- 
ton with energy E is absorbed due to K-shell ionization of 
a given element Z (rather than being absorbed by other el- 
ements or scattered by electrons) is given by the following 
expression: 



p (F \ = A Z CTK,z(£) 

Zy ' A z [a z (E) + ZcT KN (E)]+ £ a'(E) 

Z'j^Z 



(1) 



where er z - the total absorption cross section of element Z, 
<tk,z is its K-shell absorption cross section, Az - the abun- 
dance of element Z by number, ctkn - the Klein-Nishina 
cross section. For compactness of the formula we denote 
cr'(E) = A Z /[<7z'(E) + Z'ctkn(E)] - the total cross section 
due to element Z', including Thomson scattering on its elec- 
trons (see the Appendix for further details). Obviously, the 
quantity Pz(E) determines what fraction of incident pho- 
tons with energy E will contribute to the production of the 
fluorescent line of the element Z and offers a simple way to 
qualitatively investigate the dependence of the EW of the 
line on element abundances. Moreover, as it will be shown 
later, a modification of eq.JT]) gives a reasonably accurate 
method to analytically compute EWs of fluorescent lines. 

In Fig. [1] we plot the probability Pz(E) for carbon, oxy- 
gen, and iron versus energy. The left panel was computed 
for an atmosphere with solar abundances, the right panel 
for the abundances appropriate for a C/O white dwarf 
(Table [T]). As one can see from the plot, in a solar abun- 
dance case, oxygen is the dominant absorbing element for 
photon energies of up to as 7 keV. However, above the Fe 
K-shell ionization threshold of 7.11 keV, the value of Pf c 



jumps and exceeds that of Po, so that at higher energies 
the majority of photons are absorbed by iron and contribute 
to its fluorescent line. For the chemical composition of a 
C/O WD, however, the picture changes dramatically. Even 
though the K-shell photo- absorption cross-section for iron is 
larger than the oxygen cross-section at these energies, the 
increased abundance of oxygen makes it the main absorbing 
agent in the entire energy range. As a result, the Ka line of 
iron will be significantly suppressed. 

As oxygen abundance increases, Po will increase lin- 
early with Ao, insofar its contribution to the denominator 
in eq.{T} remains relatively small. However, at sufficiently 
large abundances, the oxygen term prevails and Po satu- 
rates at a value determined by the ratios of cross-sections, 
Po ~ o"K,z/(o"z-r-a"KN). On the contrary, Pf c will continue to 
decrease due to unlimited increase of a 1 in the denominator 
in eq.|T]). This behavior is illustrated in Fig. [J] where we 
plot the values of Po and Pf c estimated at their respective 
K-edges, versus Ao. As is evident from the plot, Po increases 
with Ao until the latter reaches a value of ~ 20—25 times so- 
lar value. At this abundance, oxygen completely dominates 
the opacity and nearly all incident photons that are not scat- 
tered on electrons, will be absorbed by the oxygen, regard- 
less of their energy. Further increase of oxygen abundance 
does not lead to an increase of it fluorescent line strength. 
On the other hand, the PFe curve shows unlimited decreases 
as Aq increases, asymptotically Pf c oc Aq . 



4 RESULTS 

4.1 Method of calculation 

We consider plane-parallel, semi-infinite geometry. The disk 
material is assumed to be cold and neutral. We use the 
Monte Carlo (MC) technique f o llowin g the prescription by 
iPozdnvakov. Sobol. fe Sunvaevl (|l983l ). We account for K a 
and Kg lines for all elements from Z=3 to 30 and do not 
distinguish between the K ai and K& 2 lines. Photoionization 
cross sections for given abundances are calculated using the 
analytical fits for pa rtial photoionization cross-sections from 
IVerner et alj (|l996l ). the scattering cross-section was de- 
scribed by the Klein-Nishina formula and fl uorescence yields 
for ne utral elements from Li to Zn are from lBambvnek et al.l 
(|l972h . The incident spectrum can be modeled either as a 
beam at a specific angle with respect to the normal to the 
surface or an isotropic point source above the disk surface. 
The output reflected emission can be registered at a specific 
angle or integrated over all viewing angles. In computing the 
equivalent widths of lines with respect to the total emission 



1 Note that for the purpose of this plot we fixed abundances of 
all elements except oxygen at their solar value, i.e. no condition 
of the nucleon number conservation was imposed. Such an abun- 
dance sequence does not represent any of the WD compositions 
and is employed here to investigate the behavior of equivalent 
widths in the limit of high Aq (or, cquivalcntly, high O/Fe). If 
wc use the O/Fe ratios ((O/Fe)© 26, (0/Fe) C o ~ 1.3 ■ 10 3 ) 
to characterize the oxygen overabundance, a C/O white dwarf 
would approximately correspond to the Aq/Aq,© ~ 50. 
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Figure 3. The equivalent widths of Ka lines of C, O and Fe plotted against fc/o, the fraction of H and He "converted" to C and O 
(section 14,21 cq(2)l- The left panel shows results for a power law incident spectrum with photon index of T = 1.9 and the right panel - 
for black body radiation with kT=2.5 KeV. Black solid lines show results of Monte-Carlo calculations, the dashed lines (red in the color 
version) were computed in the single scattering approximation as described in the Appendix. 



the reflected emission was mixed with the primary radiation 
as described in the Appendix. 

X-ray binaries are known to have two distinct spec- 
tral states: a bright high/soft stat e and a less luminous 
low/hard state (e.g. iGilfanovl l2010h . In the hard spectral 
state, Comptonization of soft photons on hot electrons is 
the most likely mechanism for the creation of the hard 
spectral component, which, in the energy range of inter- 
est, usually has a power law shape. In the soft state, the 
primary emission may originate in the accretion disk it- 
self or in the boundary layer on the surface of the neutron 
star ( e.g. Sunvaev fe Shakura Jl98fil : llnogamov fe Sunvaevl 
1 19991 ; iPopham fe Sunvaevl I2OO1F Correspondingly, the in- 
cident spectrum is represented by either a power law in 
which case it is characterized by a photon index T, or by 
a thermal component originating in the NS boundary layer. 
In the latter case, we approximated the boundary layer ra- 
diation with a black body spectrum with kT = 2.5 keV 
|Gilfanov. Revnivtsev. fe Molkovll2003t ). 

Along with Monte Carlo simulations, we also performed 
simple analytical calculations in the single scattering ap- 
proximation. The formulae are derived in the Appendix and 
can be used for quick analytical estimates of the strengths 
of the fluorescent lines for different abundance patterns. 



4.2 Reflection from C/O-rich material 

In order to investigate the dependence of the line strengths 
on the C/O abundance we consider a sequence of chemical 
compositions with increasing C/O abundance. Abundances 
of H and He were reduced along the sequence so that the to- 
tal number of nucleons was conserved. Abundances of other 
elements (by mass) remained fixed at their solar values, as 
well as the abundance ratio of carbon and oxygen. The po- 
sition along the sequence is defined by the parameter 



where m is concentration of the element i and mj is its 
atomic weight. Obviously, fc/o uas the meaning of the frac- 
tion of hydrogen and helium "converted" to carbon and oxy- 
gen and changes from to 1. We performed calculations of 
the reflected spectrum and computed the equivalent widths 
of fluorescent lines on a grid of values of fc/o m this range. 
The calculations were performed with Monte-Carlo simula- 
tions and using analytical expressions from the Appendix. 
The incident photons were assumed to illuminate the disk 
surface isotropically and the reflected emission was inte- 
grated over all outgoing directions. The results are shown 
in Fig. [3] for carbon, oxygen and iron K a lines, for both 
types of the primary emission spectrum. The plots show the 
equivalent widths of lines versus fc/o- 



fc/o 



n c m c + nnirir 



n H m H + n Hc m Hc + n c m c + n m 



(2) 



The overall behavior of equivalent widths is similar for 
both types of incident spectrum. It is mostly determined by 
Pz (Fig.[T]) and, as expected, the curves have shapes similar 
to those in Fig. [2] Initially, the EWs of carbon and oxygen 
lines increase nearly in direct proportion to their abundance, 
however they quickly saturate at fc/o ~ 20%. This corre- 
sponds to the oxygen overabundance of Ao/Ao,© ~ 10 in 
Fig. [2] The initial increase of the oxygen line EW is some- 
what milder than that of carbon due to screening of oxygen 
by carbon itself. In the case of the power law incident spec- 
trum, the EWs of carbon and oxygen lines remain at sub-eV 
level, even for a totally C/O-dominated composition. For a 
black body incident spectrum, however, these lines can be- 
come rather prominent at large C/O-abundances, reaching a 
value of ~ 20 eV for oxygen and ~ 2 eV for carbon. The iron 
line is the most prominent fluorescent line for solar abun- 
dances. Nevertheless, for both types of the incident spec- 
trum, its EW drops below ~ 10 eV for the C/O dominated 
composition. 
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Figure 4. The equivalent widths of Ka lines of C, O and Fe are plotted against /hoi the fraction of H "converted" to He (cf. eq[2]l. The 
left panel shows results for a power law incident spectrum with photon index of T = 1.9 and the right panel - for black body radiation 
with kT=2.5 keV. Black solid lines show results of Monte-Carlo calculations, the dashed lines (red in the color version) were computed 
in the single scattering approximation as described in the Appendix. 



4.3 Reflection from He-rich material 

We follow the same approach to investigate the dependence 
of the fluorescent lines strengths on the helium abundance 
as we used for C/0. The results are presented in Fig. 21 
where we plot equivalent widths of carbon, oxygen and iron 
lines against /He, defined similarly to /c/o in eq.([2}. 

Similar to the C/0 case, the increased abundance of 
He would have a screening effect on the elements with the 
charge Z > 2, however, the effect is much milder than in the 
case of C/O dominated material. The reason is the much 
larger solar abundance of helium and the significant drop 
of its photoabsorption cross-section at the energies corre- 
sponding to K-edges of elements which K a lines are in the 
X-ray domain. The latter factor is especially important for 
iron. Correspondingly, the overall effect is stronger for lines 
of low-Z elements, such as carbon and oxygen and is negli- 
gible for iron. For example, the equivalent widths of carbon 
and oxygen lines decrease by ~ 2 times in the He-dominated 
case. In Fig.0] one can also notice a slight, ~ 10%, increase 
of the EW of the iron relative to its solar value. This is 
a consequence of the definition of the abundance sequence, 
and does not have any physical meaning. It is caused by the 
decrease of the number of electrons, as H is being replaced 
by He along the abundance sequence. 



4.4 Realistic WD compositions 

We now compute equivalent widths of fluorescent lines 
for realistic abundance patterns expected for different 
types of white dwarfs. For these we used the results o f 
iGarcfa-Berro et all d200j ). ICil-Pons fe Garcia- Berrd (|200ll ) 
and IPandev et al.l (|2001 ) as detailed in the Section 2 and 
summarised in Table [1] The equivalent width of lines were 
computed with our Monte Carlo code assuming an isotropic 
point source above the disk surface (lamppost configura- 
tion). The reflected spectrum was integrated over all view- 



ing angles. The results of these calculations are presented in 
the Tabled 

The results are qualitatively similar to those obtained 
in the previous subsections. Perhaps a new feature is the 
appearance of the neon line in the case of the O/Ne white 
dwarf, which becomes especially prominent in the case of 
a black body incident spectrum, where its equivalent width 
can reach ~ 15 eV. 



5 DISCUSSION 

Results of the previous section suggest that absence of a 
strong K Q line of iron at 6.4 keV in the spectrum of a UCXB 
should be considered as an indication of a C/O/Ne white 
dwarfs whereas the presence of such a line points at a helium- 
rich donor. There are, however, several factors which can 
modify this simple picture. 

5.1 Gravitational sedimentation of metals 

Gravitational settling of heavy elements operates efficiently 
in the surface layers of white dwarfs, its time scale being 
considerably shorter than the evolutionary time scale of an 
isolat ed white dwarf (e.g. ISchatzmanlll958l ; IPacmette et alj 
Il98fj) . The result of the gravitational sedimentation is that 
the outer layers of the white dwarf envelope are depleted 
of metals on the time scales of ~ 1 6 yrs or shorter 
|Pacmette et al.lll98fj ; |Dupuis et ai]|l992l ). This explains the 
purity of spectra in the majority of white dwarfs, which show 
only spectral lines of hydrogen and helium. 

Extrapolating these results to the UCXB case, one may 
expect that their accretion disks should be devoid of ele- 
ments other than the main constituent of the white dwarf 
(H/He or C/O). However, although the diffusion time scale 
in the outermost layers is much shorter than the cooling 
time of an isolated white dwarf, the time scale on which the 
deeper layers of the white dwarf envelope are depleted of 
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Table 2. Equivalent width of prominent fluorescent lines for realistic white dwarf compositions. Results of 
Monte-Carlo calculations. 











Power law 








Black Body 




Element 


Energy 


Solar 


He 


C/O 


O/Ne 


Solar 


He 


C/O 


O/Ne 




keV 


















C 


0.277 


0.01 


< 0.01 


0.11 




0.35 


0.24 


8.21 







0.525 


0.25 


0.03 


0.40 


0.72 


7.10 


0.70 


13.8 


31.0 


Ne 


0.849 


0.19 


0.20 


0.06 


0.76 


3.14 


3.19 


1.00 


14.8 


Na 


1.041 


0.01 


0.01 


< 0.01 


0.23 


0.09 


0.10 


< 0.01 


3.51 


Mg 


1.254 


0.25 


0.27 


< 0.01 


0.37 


2.90 


3.17 


0.10 


4.36 


Si 


1.740 


0.84 


0.94 


0.03 


0.02 


7.23 


8.24 


0.26 


0.14 


S 


2.307 


1.13 


1.30 


0.04 


0.02 


7.19 


8.36 


0.30 


0.14 


Ar 


2.957 


0.80 


0.94 


0.04 


0.02 


3.84 


4.50 


0.17 


0.08 


Ca 


3.690 


1.03 


1.22 


0.05 


0.03 


3.60 


4.30 


0.18 


0.09 


Cr 


5.415 


0.80 


1.00 


0.05 


0.03 


1.65 


2.14 


0.11 


0.05 


Mn 


5.899 


0.52 


0.65 


0.04 


0.02 


0.98 


1.24 


0.07 


0.04 


Fe 


6.397 


76.0 


92.8 


6.58 


3.11 


127 


157 


10.9 


5.17 


Fe (K» 


7.058 


12.6 


15.7 


1.17 


0.54 


17.3 


21.6 


1.60 


0.80 


Ni 


7.470 


5.86 


7.00 


0.68 


0.33 


7.75 


9.37 


0.92 


0.44 


Ni (Kp) 


8.265 


0.93 


1.11 


0.12 


0.06 


1.04 


1.28 


0.13 


0.06 



Equivalent widths are in eV. The line energies are from lBeardenl 1119671) . The abundance patterns for different 
types of white dwarfs used in these calculations are summarized in Table [l] 



metals is rather long, much longer than the accretion time 
scale. As a result, the outer layers of the white dwarf enve- 
lope are removed faster than they are depleted of metals. 
Th e time depen dent calculations of gravitational settling 
by iDupuis et al.l show that time, required to deplete the 
outer ~ 1CP 3 (by mass) of the white dwarf exceeds ~ 10 8 
yrs. On the other hand, for the accretion rate of ~ 10 -9 
Mo/yr, the outer layer of mass AM ~ 10 -3 Mq is removed 
within At ~ 10 6 yrs. Therefore we do not expect that grav- 
itational settling of heavy elements is an important factor 
determining the composition of accretion disks in UCXBs. 
This conclusion is confirmed by detection of iron lines in 
the spectra of UCX Bs with He-rich donors (|Asai et al 1 l200d ; 
iBoirin et ai]l2004 ). 



5.2 Ionization state of the accretion disk 

The material on the surface of the accretion disk may be 
ionized due to internal heating as well as due to irradiation. 
This can significantly modify the shape of the reflected emis- 
sion and strengths of fluorescent lines. For example, if carbon 
and oxygen are fully ionized, they do not contribute to the 
ionization cross-section and the fluorescent line of iron will 
have its nominal, ~solar abundance strength even for the 
C/O-dominated disk. On the other hand, partial ionization 
of oxygen would lead to appearance of much stronger lines 
of OVII-OVTII, but would not eliminate the effect of screen- 
ing of elements with higher charge Z. Correspondingly, the 
iron line would remain suppressed in this case. 

Illumination of the accretion disk by X-ray photons 
produced in the vicinity of the compact object (e.g. in 
the boundary layer) results in the appearance of the ion- 
ized skin - a thin surface layer of hig hly ionized material 
|Navakshin. Kazanas fe K allrnan 2000). Beneath this layer, 



the disk remains in the ionization state determined by heat- 
ing due to viscous dissipation. The Thomson optical depth r 
of the ionized skin is determined by the "gravity parameter" , 

H L^dd 



A = 



R I/ x cos ( 



(3) 



where H/R is the aspect ratio of the Shakura Sunyaev ac- 
cretion disk, Lx is the luminosity of the irradiating source, 
^Edd - Eddington luminosity for a neutron star and 9 - 
the angle of the incident radiation with respect to the nor- 
mal to the disk surface. This parameter characterizes the 
strength of gravity relative to that of radiation pressure. 
For typical parameters of UCXBs, H/R ~ few x 10 -2 , 
L^/L-Edd ~ 10~ 2 — 10 _1 a nd assuming cos 6 ~ 0.1 we obtain 
A ~ 3 — 30. According to iNavakshin. Kazanas fe Kallmanl . 
this value corresponds to low to moderate illumination, in 
which case the Thomson optical depth of the ionized skin 
is small and it does not significantly distort the spectrum 
of radiation reflected off the inner layers of the disk. The 
ionization state of these layers is determined by the viscous 
dissipation. 

The effective temperat ure of the Shakura-Sunya ev disk 
due to viscous dissipation jShakura fe Sunvae v 1973) is 



Teff = 



/ 3GMM \ 1/4 /r, 
V SWo 3 J 



3/4 



1/4 



(4) 



where a = 5.67 ■ 10 _5 erg cm -2 s _1 K -4 is the Stefan- 
Boltzmann constant, M is the mass of the accretor (~ 
1.4M for a neutron star), M the mass accretion rate and 
r is the inner radius of the disk, which in the case of the soft 
state of a neutron star is equal to its radius. For the mass 
accretion rate in the M ~ 10~ 10 — 10~ 8 M /yr range (cor- 
responding to L x ~ 10 36 1 n3S 
temperature ranges from 



10 38 ergs 1 ) the maximum disk 
2 ■ 10 6 to ~ 8 ■ 10 6 K. Obviously, 
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hydrogen and helium are fully ionized throughout the most 
of the accretion disk at any value of the accretion rate rel- 
evant to UCXBs (as they should be, in order for accretion 
to proceed in the high viscosity regime). Iron is virtually 
never fully ionized, its ionization state being below FeXXIII 
in the entire temperature range of interest, T c g < 8 • 10 6 
K. The ionization states of carbon, oxygen and neon, how- 
ever, depend critically on the mass accretion rate. Indeed, 
for plasma in the collisional ionization equilibrium, in coro- 
nal approximation, 90% of carbon (oxygen) is in the fully 
ionized state at temperatures of ~ 1.6 ■ 10 6 (~ 3.5 • 10 6 ) K 
(|Shull fc van Steenberdll982T ). Therefore, even at the mod- 
erate luminosity of ~ 10 37 erg/s, carbon is expected to be 
ionized in a significant part of the inner disk, out to r ~ 10r o . 
Oxygen and neon, on the other hand, are fully ionized in 
the inner disk at high luminosities, ~ 10 38 erg/s, but re- 
main only partially ionized at moderate and low luminosi- 
ties, L x ~ 10 36 - 10 37 erg/s. 

Thus, one should expect luminosity dependence of the 
iron line strength in UCXBs. At low and moderate lumi- 
nosities, it is determined by the chemical composition of the 
accreting material, as described in the previous section. In 
particular, it remains at the solar value in systems with He- 
rich donors and is suppressed in the case of a C/O/Ne donor. 
As luminositj0 approaches ~ 10 38 erg/s, the iron line should 
recover its nominal strength, determined by the abundance 
of iron in the accreting material. 



5.3 Reflection from the white dwarf surface 

For a binary system consisting of a 1.4M© NS accretor and 
a O.7M WD donor, the Roche lobe of the WD subtends 
a solid angle of ~ 0.3 sr as viewed from the neutron star. 
In the idealized lamppost geometry the solid angle of the 
accretion disk is ~ 2n. Therefore, the donor star does 
not contribute significantly to the total reflected emission. 
However, the lamppost geometry is obviously a too crude 
approximation. In more realistic geometries, taking into ac- 
count the geometry and emission diagram of the primary 
emission, the role of reflection from the donor star may no 
longer be negligible. 

For a T ff ~ 10 4 K white dwarf, helium is expected to 
be partly or fully ionized, wher eas carbon, oxygen and iron 
should be only partly ionized (|Dupuis et al.lll992f ). There- 
fore, for a He white dwarf, the fluorescent line of iron should 
have its nominal strength. For a C/O white dwarf the iron 
line should be suppressed due to screening by carbon and 
oxygen, as described above. Depending on the exact value 
of the temperature in the white dwarf photosphere, carbon 
and oxygen may be in the medium to high ionization state. 
In the latter case, one may expect appearance of strong res- 
onant lines of their highly ionized species. Among these, of 
special interest would be lines of oxygen. 



2 Determination of the exact value of luminosity at which the 
transition happens, requires detailed calculations of the ionization 
state of the material at the surface of the accretion disk and is 
beyond the scope of this paper. 



5.4 Iron abundance 

In order to determine the iron abundance, we relied on the 
assumption that the total mass of elements, heavier than 
O/Ne, in the donor star, does not change. This was moti- 
vated by the fact that nuclear reaction chains do not involve 
heavy elements. Moreover, the heavier elements were fixed at 
their solar mass fractions, in order to make comparisons be- 
tween spectra originating in H-rich and H-poor atmospheres. 
These assumptions, along with the requirement to conserve 
the total number of nucleons, effectively limited the maximal 
C,0,Ne/Fe abundance ratios. This in turn limited the mini- 
mal value of the iron line EW at small but moderate values 
of ~ several eV (Fig. Table [21 . However, if one allowed 
further increase of the O/Fe abundance ratio, the equiva- 
lent width of the iron line can decrease to arbitrarily small 
numbers, in inverse proportion to this ratio, as illustrated 
by the Fig. [2] 



5.5 Comparison with observations 

Although we plan to perform detailed comparison of our 
predictions with observations of UCXBs with Chandra and 
XMM-Newton in a follow-up publication, we present some 
preliminary conclusions based on the published results. As 
already mentioned in the introduction, an iron li ne has been 
detected in the spectrum of UCXB 4U 1916-05 (|Asai et all 
|2000| ; iBoirin et all 12004!') . which is believed t o have a He- 
rich donor ( Nelemans. Jonker fc Steeghsl2p06f). Simil a rly 4U 
1820-30, for which iBildstenl I | l995l ) and ICummind (|2003T l 
have predicted a He-rich do nor, has also been f itted with 
an iron K a emission line (jCackett et all l201Ch . On the 
other hand, objects 4U 0614+091 and 4U 1543-624, that 
show si gnatures of C/O-rich material at optical wave- 
lengt hs (|Nelemans et alj|2004l ; [Nelemans. Jonker fc Steeghsl 
|2006[). seem to have very faint, if any, iron K a emission 
(|Madei et alll2010l ; iMadei fc Jonkerll201ll ). 

However, the picture is far from being complete and 
conclusive, with different analysis and different observa- 
tions givi n g som etimes contradicting results. For example, 
iNg et all (|201Ch found an iron line in the XMM spec- 
trum of 4U 1543-624 with an EW of ~ 30 eV, whereas 
ICostantini et all (|2012T ) did not detect an i ron line in the 
of XM M spectra of 4U 1820-30. In addition iKuulkers et all 
interpret the observed bursting activity of 4U 
0614+091, by considering a He-rich or hybrid donor, instead 
of a C/O-rich one. We plan to conduct systematic investiga- 
tion of XMM-Newton and Chandra observations of UCXBs 
in a separate paper. 



6 SUMMARY AND CONCLUSIONS 

We have shown that non-solar composition of the donor star 
in ultra-compact X-ray binaries may have a dramatic effect 
on the reflected spectral component in UCXBs, significantly 
modifying the strength of fluorescent lines of various ele- 
ments. 

To identify the main trends, we considered an idealized 
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case of an optically thick slab of neutral material with signif- 
icantly non-solar abundances. We considered two abundance 
patterns, corresponding to a He and C/O white dwarf. We 
found that from the observational point of view, the non- 
solar composition of the reprocessing material most pro- 
nouncedly affects the strength of the fluorescent line of iron 
at 6.4 keV. Although increase of the carbon and oxygen 
abundances does lead to some increase of the strengths of 
corresponding fluorescent lines, their equivalent widths sat- 
urate at the (sub-)eV level due to the effect of self-screening. 
On the other hand, the equivalent width of iron decreases 
nearly in inverse proportion to the C/O abundance and the 
line is expected to be significantly fainter for the chemical 
composition of the C/O white dwarf. This is caused by the 
screening of the presence of heavy elements by oxygen. In 
C/O-dominated material, the dominant interaction process 
for a E > 7 keV photon is absorption by oxygen rather than 
by iron, contrary to the case of solar composition. Screening 
by helium is significantly less important, due to its lower ion- 
ization threshold. Moreover, helium is expected to be fully 
ionized in the accretion disks of UCXBs. Consequently, in 
the case of He-rich reprocessing material, fluorescent lines 
of major elements are near their nominal, solar abundance 
strength. Thus, the equivalent width of the fluorescent line of 
iron can be used for diagnostics of the donor star in UCXBs 
by means of X-ray spectroscopy. 

In the realistic case of reflection in UCXBs, gravita- 
tional settling of elements and ionization of the disk ma- 
terial may, potentially, complicate and modify this picture. 
Simple comparison of the diffusion time in white dwarf enve- 
lope and the accretion time scale suggests that gravitational 
settling is not fast enough to deplete iron and other heavy el- 
ements in the accreting material. However, a more accurate 
consideration of the physical state of a Roche-lobe filling, 
mass losing white dwarf may still be needed for the final 
conclusion. 

Ionization of the disk material at high mass accretion 
rate may lead to luminosity dependence of the discussed 
effects. In particular, as oxygen in the inner parts of the 
C/O-dominated disk becomes fully ionized at high mass ac- 
cretion rate, its screening effect vanishes and the iron line 
should be restored to its nominal value, determined by the 
abundance of iron in the accretion disk. Comparison of ion- 
ization curves with the effective temperature distribution in 
the Shakura-Sunyaev accretion disk suggests that it should 
happen at the log(Lx) ~ 37.5 — 38 level. At lower luminosi- 
ties, log(Lx) < 37, oxygen in the inner disk is in the low- to 
medium ionization state and the idealized picture outlined 
above holds, at least qualitatively, and the strength of the 
6.4 keV iron line may be used for diagnostics of the nature of 
the donor star. In particular, its absence points at the C/O 
or O/Ne white dwarf, while its presence suggests a He-rich 
donor. 
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APPENDIX A: SINGLE SCATTERING 
APPROXIMATION 

We consider reflection from a semi-infinite atmosphere in the 
single scattering approximation. The spectral intensity of re- 
flected emission S TC n (phot sec -1 cm -2 keV" 1 sr" 1 ) is given 
by the following expression. 



oo 

J dz SeC 6» in g-" CT tot(E)zsec£) in [E)P BC (ftm, flout ) 



e -"Otot(-E)zscce out (Al) 

where Sp r (E, 2) is the spectral intensity of the primary ra- 
diation, axis z is normal to the surface of the atmosphere 
and is directed inwards, 8Q ln is an infinitesimal solid angle 
around the direction of incidence ni n , 9 m is its polar an- 
gle with respect to the axis z. The spectral intensity of the 
reflected emission is computed at the direction n ou t, which 
polar angle is out . P S c(ni n ,n ou t) is the probability that the 
photon, entering the medium from the direction ni n is scat- 
tered in the direction n out , it is normalized so that: 



-fsc\Hin, flout O 



(A2) 



The n is the density of the material and CTtot = fabs + o" sc is 
the total cross section. The absorption cross section due to 
photoionization <r a bs is given by the following expression 

30 

<Tab s (£) = A z az(E) (A3) 

Z=l 

where we account for all elements from Z = 1 to Z = 30, 
Az is the abundance of element Z by the particle number, 
and (Jz is the photoionization cross-section for all shells of 
element Z. It is cal culated using the second version of the 
IVerner et all (|l996h subroutine. (j sc is the scattering cross 
section per hydrogen atom, given by 



a sc (E) = J2 ZA ^t 



(A4) 



where ctt is the Thomson cross section. Note that we con- 
sider Compton scattering in the low energy limit and ignore 
change of the photon frequency during scattering in deriving 
eq. (|AT) . 

For a semi-infinite atmosphere, the reflected spectrum 
does not depend on the density of the material, only on its 
chemical composition. 

We assume for simplicity that the energy and angular 
dependencies of the primary radiation can be factorized: 



S pr (E,n) = S {E)P pl (n) 



(A5) 



where the P pr (ni n ) describes the angular distribution of the 
primary radiation and is normalized so that 



Ppr(n) dtt : 



(A6) 



SS IE a(E, n ut) = S pv (E, ni n ) 8£li n 



and So(E) is proportional to the total luminosity of the 
primary emission and has units of phot sec -1 cm" 2 keV -1 . 
We evaluate the reflected emission within solid angle Af2 ut 
around direction of interest n ou t: 

Srea(E)= I S Icfi (E,n out )dn (A7) 

Integrating eq. Al over ingoing and outgoing directions and 
from z=0 to z = oo we obtain the following expression for 
the reflected continuum S Ta a 

a BC (E) 



S Ic n(E) — So(E) 



a ac (E) + o-abs(-B) 



(A8) 
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where the factor R Te a accounts for the geometry of the prob- 
lem and is given by 

-Rrefl = / dQi n P pl (n in ) / df2 out 

sec 0i n (secftn + sec6< ou t) _1 P S c(ni n , n out ), (A9) 

For Thomson scattering P sc (ni n , n out ) depends only on the 
scattering angle 9 SC and is given by the standard Rayleigh 
formula 



sc \ w sc 



) = 



3 (1 + cos 2 



(A10) 



(All) 



8tt 2 
Ignoring angular dependences, 

For the case of isotropic incident radiation and reflected 
emission integrated over all outgoing angles (the geometry 
assumed in Fig. [3l [4j) eq. A9 yields R ro fl = 1/8. 

Using the same approach and taking into account the 
fluorescent yield, we can compute the fluorescent line flux 
Piinc (phot sec" 1 cm -2 sr): 

oo 

5Fi ine (n out ) = <5f2i n j dE S P r(E,n in ) 



OC 

/ 



dz sec Si n e 



-n<Ttot (E)zsec 



A z n aK,zYzPiine(n ou t) 



-"CT tot (B lins ,)2 sec 6 out 



(A12) 



In the above expression, ok.z, is the K-shell absorption cross 
section of element Z and Yz is the fluorescence yield of its 
K a line, -Eii ne is the K a line energy and Ek is the energy 
of the K-edge. Pnne(n ut) is the angular distribution of the 
fluorescent emission, assumed to be isotropic (Pii n c(n ut) = 
l/47r). Integrating over all angles and over z — — > oo we 
obtain 



+ [o"tot(-B) 2 - (Ttot(i?line) 2 ] ln[tT t ot(£ ; ) + 0tot (Pline)] 



and G2 is given by, 

G2 = 8(Jtot(P) 2 Crtot(Plinc) 



-a tot (EYln[a tot (E)] (A16) 



(A17) 



To characterize the strength of emission lines we eval- 
uate their equivalent widths with respect to the total con- 
tinuum emitted within the solid angle Af2 ut around the 
direction of interest. 



EW 



Fun 



(A18) 



Stot (-En 

no I 

The total continuum includes both the reflected continuum, 
given by eq. (|A8[l and the fraction of the primary continuum 
emitted in the solid angle Afi ou t 



Stot{E) = S re a(E) + S pr (E) 
where, similarly to S' ro fl(E), 



S pr (E) = 



Spr (E,n out ) dQ out 



(A19) 



(A20) 



Thus, using equations I A8I and I A9I for the reflected con- 
tinuum, eqs. IA13I and IA14I for the fluorescent line flux, the 
equivalent width of the fluorescent line can be computed 
from eqs. IA181IA20I Comparison with Monte-Carlo calcula- 
tions show that the single scattering approximation works 
nearly perfectly at low energies, E < 2 keV, where absorp- 
tion dominates scattering. At higher energies, multiple scat- 
tering becomes more important, resulting in a ~ 10% offset 
between analytical and Monte-Carlo results for the iron K a 
line. This is further illustrated by Figs. [3] and [4] showing 
dependence of equivalent widths of various fluorescent lines 
on the chemical abundances computed in single scattering 
approximation and using the Monte-Carlo code. 



Fii„c = J S (E)G(E)dE 

Ek 



(A13) 



where as before, tilde denotes integration over solid angle 
ASlout and G(E) is given by the following integral 



G(E) 



dfli n P pr (8in) I dfiout-Pline(Wout) 

ASl in JAn out 

Y z Az a K ,z(E) sec in 



, crtot (E) sec din + ftot (Ehob) sec out J A14 ^ 

A similar approach was used by IChurazov et al.l (|2008h in 
calculating the Earth albedo. Similarly to R Ten , for the 
geometry of Fig. [3l [4] (semi-infinite slab illuminated by 
isotropic incident radiation, reflected emission integrated 
over all outgoing ang les), eq. lAlU can be integrated ana- 
lytically to give: 

G(E) = YzAz a K , z (E) — (A15) 
Where Gl is given by, 

Gl = CTtot(-Elinc)[o-tot(-E') + Cr t ot(Pline)ln[(Jtot(-Elino)]] 
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